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Abstract 

A Massive Compact Halo Object (Macho) that lenses a background star will magnify 
different parts of the rotating stellar disk by varying amounts. The differential magni- 
fication will cause a shift in the centroid of the star's spectral lines during the lensing 
event. The shift is proportional to the ratio of the stellar radius to the projected 
separation of the Macho from the star. It therefore provides a direct measure of the 
Einstein ring radius, and so also a measure of the Macho's proper motion (angular 
speed). This measurement can remove some of the degeneracy between mass, distance 
to the lens, and transverse velocity that exists in the interpretation of results from 
ongoing microlensing experiments, and is an independent test of the lensing nature of 
the event. We show that using the high precision attainable by stellar radial velocity 
measurements, it is possible to measure proper motions for ~ 10% of Machos that 
lens A-stars in the Large Magellanic Cloud (LMC), i.e. ~ 7% of the type of relatively 
high-magnification events that have been reported to date. If this proper-motion 
measurement were combined with a parallax measurement of the "reduced velocity" , 
then the Macho mass, distance, speed, and direction could each be separately de- 
termined. The shift can be measured for ~ 20% of the A-star events generated by 
Machos in the dark halo of the LMC. This in turn would provide a measurement of 
the fraction of LMC vs. Galactic Macho events. 

Subject Headings: spectroscopy - dark matter - gravitational lensing - Magellanic 
Clouds 



1 Introduction 

Two independent groups have recently found a total of three candidate microlensing 
events, apparently caused by Massive Compact Halo Objects (Machos ) along the line 
of sight toward the Large Magellanic Cloud (LMC) (Alcock et al. 1993; Aubourg et 
al. 1993). The events are achromatic, have maximum magnifications of Amax = 6.8, 
3.3, and 2.5, and characteristic times oo^^ = 17 days, 30 days and 26 days. The light 
curves fit the form first predicted by Paczyhski (1986) quite well: 

A[xit)] = ^^^^|t7^, xit) = ^^- + u\t-t,Y. (1) 

Here to is the midpoint of the event, uo~^ is the characteristic time, and j3 is the 
dimensionless impact parameter (normalized to the Einstein ring radius). 

Of the three measurable parameters in equation (1), only the characteristic time, 
oo~^ yields any information about the Machos. It is related to the underlying physical 
parameters by 

_i _ e^DoL _ [4GMDol(1 - /^ol//^lmc)]'/' ,^, 

V vc 

where 6^, is the Einstein ring radius, M is the Macho mass, v is its transverse velocity, 
and Dql is its distance from the observer. The distance of the lensed star is denoted 
-Dlmc- For any given event, one cannot determine the mass, distance, and speed 
separately. For an ensemble of events, the typical speed and distance are expected 
to be f ~ 200 km s^"*^ and Dql ~ 10 kpc, respectively. Hence, M ~ (7Oc<jdays)^^M0. 
If the first few events prove typical, it will imply a mass scale of ~ 0.1 Mq, near the 
hydrogen-burning limit. 

One would like to determine as much as possible about the distribution of Machos, 
not just the mass scale. For example, one would like to know if the Machos actually 
lie in a halo as opposed to a disk or thick disk, whether the halo is spherical or 



oblate, whether the halo is truncated or extends all the way to the LMC, whether it 
is rotating, whether Machos populate the LMC as well as the Galaxy, and what the 
distribution of Macho masses is. To extract these additional pieces of information, 
new methods of analyzing the data as well as new experiments are required. 

Gould (1992, 1993, 1994a,b), Sackett & Gould (1993), and Gould, Miralda-Escude 
& Bahcall (1994) have developed a number of different techniques for extracting 
additional parameters from lensing events. In particular, Gould (1994a) showed that 
the non-zero size of the lensed stellar disk modifies the light curve of the event in a 
manner that allows one to measure the Macho's Einstein radius, ^*, and hence its 
proper motion (angular speed), uO,,. In practice, however, the effect is measurable 
only for a fraction of events during which the Macho passes over the line of sight to 
the face of the source star, (1/550)(M/0.1 Mo)^^/^ for Machos in the Galactic halo 
and (l/7O)(M/O.lM0)-i/2 ^^ Machos in the LMC. 

In this Letter, we predict a second, spectroscopic, effect which also arises from the 
non-zero size of the source, combined with the source rotation. Differential magnifica- 
tion of the stellar disk during the event induces a shift in the stellar spectral lines. By 
measuring the shift as a function of time one can determine the Einstein ring radius 
and so the proper motion of the Macho. In contrast to the photometric effect analyzed 
by Gould (1994a), this spectroscopic effect can be measured even when the Macho is 
many stellar radii from the source. Hence, the proper-motion measurement can be 
made for a larger fraction of events. We estimate that with an 8m telescope the effect 
can be measured in ~ 7% of photometrically detected relatively high-magnification 
{j3 <) 0.5) Macho events. The effect can be measured in ~ 15% of high-magnification 
events generated by Machos in the LMC. 

Measurement of the proper motion of Galactic Machos would, in itself, remove 
some of the degeneracy among the three Macho parameters, mass, distance, and 



transverse speed. However, if this measurement were combined with a parallax mea- 
surement of the "reduced velocity" , then the degeneracy could be completely broken 
and the three parameters plus the transverse direction could be separately determined. 
The proper motion of Galactic Machos is ~ 15 times greater than that of LMC 
Machos. Hence, a proper-motion measurement clearly distinguishes between the two. 
By identifying even a handful of LMC Machos, one could determine the fraction 
of events generated by them. This would both provide direct information about the 
LMC halo and remove a major background to the primary Galactic signal (see Sackett 
& Gould 1993). 

2 Spectral Shift in Lines of a Microlensed Star 

The spectral line profile of an unlensed star is broadened by the effect of the star's 
rotation. Different parts of a finite stellar disk are magnified by varying amounts 
during a Macho lensing event. In order to determine the time- dependent rotational 
line profile, we apply equation (1) to calculate the magnification at each point on the 
star. Expressing all angles in units of the Einstein radius, ^*, the line profile is 



f{v,t) = C r'^dz2A{\^{t) + z|)5(z), V = ^t;rot + vo (3) 

where p is the stellar radius, x(t) is the projected separation of the stellar center from 
the Macho, 5'(z) is the surface brightness as a function of position, z, on the stellar 
disk, fj-ot is the projected rotation speed of the star, Vq is the velocity of the stellar 
center of mass, and C is a constant. The Z2 axis is defined to be aligned with the 
stellar rotation axis. The effect of limb-darkening can be approximated (e.g. Allen 
1973) by the form: 

5(z) = 1-U2-V2 + U2il - |z| Vp2)i/' + t;2(l - |z| Vp'), (4) 

with U2 = 0.99 and V2 = —0.17 at 4500A. We assume the star is devoid of star-spots. 
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The magnitude of the separation function x{t) = J (3"^ + (cut)^ is known from the 
overall light curve. One can integrate equation (3) numerically and compare the 
expected profiles with the rotational profiles observed in a series of measurements. It 
will generally be more practical to measure the shift in the centroid of the line, i.e. 
the star's apparent radial velocity, as a function of time: 

^ Sdvf{v,t)v 
Jdvf{v,t) ■ 

Apart from the rotational broadening described by equation (3), stellar line profiles 
can be intrinsically broadened by thermal and turbulent motions, by the natural line 
width, and by pressure broadening. The intrinsic broadening does not, however, affect 
the line centroid. 

To study the behavior of the velocity shift, we first make the simplifying assump- 
tions that p '^ X '^ 1 and S{z) = constant. The integral for the line profile (eq. 3) 
can then be approximated as 



"^rot 

and the line centroid (eq. 5) becomes 



V — vq P cos a + ut sin a 



1 + -p 



V,ot P^ + itOt) 



(6) 



^ pcosa + utsina 

{vit))=vo + -v...p ^,^^^^^, , (7) 

where ^ = 1 and a is the angle between the Macho velocity d'^/dt and the star's 
projected rotation axis. For a limb-darkened stellar profile as described by equation 
(4), equation (7) remains valid with S, = (30 — 14-U2 — 20f2)/(30 — 10^2 — 15^2), e.g., 
^ = 0.86 at 4500A. We find numerically that for p <, x <^ 1, the velocity shift in 
equation (7) should be reduced by a factor [1 — 0.2(p/a;)^]. 

Figure 1 shows the shift in the line centroid, A{v{t)) = {v{t)) — fg as a function of 
time for an event with timescale uj~^ = 17 days and a star with fi.ot=100 km s^^. The 
dependence of the line-shift curve on the parameters p, f3, and a, the angle between 
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the stellar axis and the Macho trajectory, is illustrated in the separate panels. As 
reference, for a O.IMq Macho in the Galactic halo at a distance of 10 kpc lensing 
an A star of radius 2.5-Rq in the LMC, p = 0.001. Such a Macho in the LMC halo, 
lensing an LMC star 2.5 kpc behind it, will have p = 0.01. The fraction of lensing 
events having dimensionless impact parameter of (3 or less (and hence maximum 
magnification /5~^ or greater) is simply /5. All values of the angle a are equally likely. 

As can be seen from equation (7) and Figure 1, varying the value of p, the ratio 
of the apparent stellar radius to the Einstein radius, affects only the amplitude of 
the line-shift curve. Changes in the impact parameter f3 affect both the amplitude of 
the curve and its shift in time relative to the time of maximum total magnification. 
Changes in the projected angle a between the stellar spin-axis and the Macho tra- 
jectory affect the amplitude, the time-shift, and the shape of the curve. The impact 
parameter f3 can be fit independently to the photometric light curve. The parameters 
p and a can then be fit unambiguously to an event, given enough signal to noise, and 
the one interesting parameter, p, measured. Since the physical stellar radius is known 
from the star's spectral type, as is the distance to the LMC, p provides a measure of 
the Macho Einstein ring radius, 6^ and the Macho proper motion uj9^. Knowledge of 
the proper motion removes one part of the degeneracy between Macho mass, distance, 
and transverse velocity in the interpretation of the photometric light curve (see eq. 2). 
For example, for an assumed velocity v, the distance Dql to the lens can be derived, 
and then from this distance and 6^, the mass M can be calculated. 

From equation 7 and the relevant range of parameters, a velocity resolution of 
(5f /frot ~ p ~ 0.01-0.001 is required in order to measure p in some fraction of lensing 
events. A temporal resolution of order hours is needed. We examine the observational 
aspects more carefully below. 

The stars being monitored by the MACHO and EROS collaborations are ~ 1/3 



main-sequence stars, mostly A stars, and ~ 2/3 red giants (C Alcock 1993, private 
communication). The projected rotation speeds of A stars have a broad distribution 
with a mean fi.ot=150 km s"-*^ for AO types, decreasing to frot=80 km s^-*^ for FO (e.g. 
Allen 1973). Giants bluer than about G2 have similar rotation properties, with mean 
"^rot ~ 75 - 55 km s~^ for FO III to GO III. A "rotation boundary" exists for giants at 
spectral type G3 (Gray 1989; a similar boundary exists in the main sequence at F5). 
Giants redder than G3 rotate very slowly, with mean velocities of f rot ~ 5 - 2 km s~^ 
for G3 III to K3 III. Giants later than G3 rarely or never rotate faster than 9 km s^^. 
Furthermore, their line broadening is dominated by photospheric macroturbulence 
of order ~ 4 - 8 km s^^ (Gray 1989), which affects the accuracy with which their 
frot can be determined. (Below, in § 3, we will show that the accuracy of the proper- 
motion measurement is fundamentally limited by how well fj-ot is known.) A precision 
5v ~ 100 m s~^ is required for A stars and ~ 5 m s~^ for red giants. 

Such precisions can be obtained today in stellar radial- velocity measurements. For 
example, the stellar radial- velocity survey at the CfA achieves accuracies ~ 1 km s^^ 
for 12 mag A stars and ~ 200 m s^^ for G-stars in 20- minute exposures on a 1.5m 
telescope (Latham 1992; T. Mazeh 1993, private communication). The survey uses an 
intensified reticon detector, which can utilize only one echelle order. Since the survey 
is not optimized for A stars, the one order includes only weak metal lines, and none 
of the Balmer lines which are strong in A stars. Modern CCD detectors can have 
at least 5 times higher quantum efficiency. CCDs can simultaneously record many 
echelle orders and so increase the observed number of lines and hence the precision. 
The precision of radial velocity measurements increases as (No.photonsxNo.lines)^/^, 
but is ultimately limited by systematic effects, such as spectrograph stability. It 
is plausible that a 17 mag A star in the LMC could be measured to ~ 100 m s^^ 
in a one-hour exposure with a CCD on an 8m telescope with a spectroscopic setup 



optimized for A stars. (Recall that the stars being monitored, with typical unlensed 
magnitudes of 19-20, are in the process of being magnified by several magnitudes.) 
While improved velocity resolution is possible for red giants because of their numerous 
and strong lines, techniques yielding resolutions of < 10 m s^^ have been applied only 
to very bright (V ~ 5 mag) stars (e.g. Walker 1992), and may not be feasible here. 
Because of this difficulty as well as the previously mentioned problem of determining 
the rotation speed of red giants, it may not be possible to make a spectroscopic 
determination of the proper motion of Machos which lens this class of star. When we 
estimate the fraction of measurable events below, we therefore consider only A stars. 

3 Error Estimates 

We examine now the errors that would be involved in an actual measurement. Pa- 
rameterizing the separation x by an angular variable 6, 

X = Psec9, (8) 

equation (7) becomes 

c c ^ 

(v) = Vq -\ — ^^pcosctcos^^ H — ^^p sin q; COS ^ sin ^ = ^^aifiiO). (9) 

where 

ao = — , ai = psmQ;, a2 = pcosQ;, (10) 

h{e) = Q, fi{e) = Qsmecose, h{e) = Qcos'9, (ii) 

and Q = frot'C/4/9- (Note that by including Q in the trial functions, fi{6), we have 
implicitly assumed that Q is known. If there is a fractional uncertainty in Q - due 
for example to an uncertainty in f rot ^ this will induce a similar uncertainty in p over 
and above those analyzed below.) 
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We now suppose that a series of measurements are made of the velocity centroid 
and that the errors in these measurements scale inversely as the square root of the 
product of the exposure time and the luminosity. We normalize these to the error 
of a one hour observation at maximum magnification, which we take to be (Tq = 6v. 
Hence for a general observation of length 6t and position 6, 

a\e) = (Sv?—sece. (12) 

ot 

We note that during the time interval dt, the lens moves by d6 = cos^ 6uj/(3dt, where 

cu^^ is the time required to cross an Einstein ring radius. According to standard linear 

theory (Press et al. 1989), the inverse covariance matrix bij for the three parameters 

tti is given by 

For purposes of illustration, we assume that the observations are conducted continu- 
ously between ^min and ^max- Then, the above sum can be converted into an integral 

n ( sec Q sin Q cos Q \ 

bij = A'^ rf^ sin^ sin^^cos^ sin^cos^^ , (14) 



where 



cos 9 sin 9 cos^ 9 cos' 



4:6V 



3, 



A = -Jf3ujhi. (15) 

For example, for an event with peak magnification of 5, a; = (17days)^^, and 8 hours 



observations per day, A = 45f/(frotO\/0-2(17)~^/8. 

For simplicity, we consider the case where the observations are made symmetrically 
around the peak of the light curve, 6'max = — ^min = ^o- The total length of the 
observations is then T = 2{j3/uj) tan^o- For this case, the even and odd components 
decouple and bij takes a relatively simple form. We invert b to find the covariance 
matrix c = b^^. The relevant terms are 

Cii = - CSC^ ^0 A^ C22 = - sm^o ^- . ^ ,. ,^ r-- A^, (16) 



2 ^ ' ^^ 2V 3 In tan(^o/2 + 7r/4) 
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for the variances of ai = psina and 02 = pcosa, respectively. Note that the quantity 



of interest is p = Ja\ + a^- 

In the above formal treatment, we have implicitly assumed that the source star is 
moving with constant velocity relative to the Earth. However, in general there will 
be some relative acceleration due to the motion of the Earth about the Sun and to 
the possible motion of the source about an unseen companion. The Earth's motion is 
known and can be taken out. The motion of the source could likewise be taken out if 
it were measured by long-term observation of the source radial velocity. However, it 
is much easier to measure the source velocity when the source is highly magnified, so 
the best constraint on the acceleration of the source may come from the observations 
during the lensing event. This acceleration will be very nearly linear provided that 
the orbital period is much larger than the span of observations. We assume that this 
is the case. 

The effect of an acceleration, g, is to add a term which is linear in time, i.e. 
oc tan^. Formally we write f^{0) = QtanO; 03 = g/Q. Once this additional degree 
of freedom is allowed, the covariance matrix becomes four-dimensional. However, 
under the assumption of a symmetric interval of observations, the even and odd 
components decouple. We can therefore examine two smaller submatrices, b'^ for 
components and 2, and b^ for components 1 and 3. For the odd terms, we find 

, _ _ ._2 f^" ,^ f siv? OcosO sin ^ tan ^ \ , . 

^'^-^ J.eo [smOtane tan^^sec^j' ^^^' 

which implies a variance of ai = p sin a of 

^_ ^ 1 / sin-^gp _ ^ [lntan(go/2 + 7r/4)-singo]^ \ -^ ^^ 



^^^ 2V 3 "tan^osec^o-ln tan(^o/2 + 7r/4). 

The variance C22 of 02 = pcosa remains as in equation (16). The rms uncertainty in 



p (averaged over all angles a) is 6p = [(cn + C22)/2] 



11 



1/2 



Figure 2 shows 5p/ p as a function of the total time span of the observations T, 
where we have assumed u'^ = 17 days and Sv/vj-otC, = 10^^ per one-hour observation 
at maximum magnification. The two panels show the errors in p for p = 0.001 (typical 
Galactic Machos) and for p = 0.01 (typical Machos in the LMC) for various values 
of the impact parameter j3. We have assumed 8 hours of observation per day in the 
Galactic Macho panel (p = 0.001), and just one hour per day for the LMC halo 
Machos (or, equivalently, a longer observation on a smaller telescope, achieving the 
same Sv/vj-ot)- Assuming one is interested in < 30% error in p (i.e., a 3o" detection) 
and one can dedicate < 8 days of telescope time to the monitoring of the Macho 
event, we see that /3 < 0.1 is required for typical Galactic events. One measurement 
per night for 8 days would be sufficient to detect LMC Machos with (3 < 0.2. If the 
LMC has a dark halo, ~ 10% of all events are expected to be LMC events (Gould 
1993). Hence ~ 10% of all A-star events can be detected at the 3 a level. Note that 
the detection capabilities of the MACHO and EROS experiments are biased toward 
low-/3 events, and that all their detected events so far have magnification > 2 (i.e. 
f3 < 0.5). The velocity shift is therefore measurable in ~ 20% of the A-star events of 
the type that is being detected, or ~ 7% of the events being detected on all stars. 

4 Discussion and Conclusions 

We predict the existence of a spectroscopic line-shift, resulting from differential mag- 
nification of a rotating stellar disk, during microlensing events in the ongoing MACHO 
and EROS experiments. The effect can be measured and monitored on large tele- 
scopes for a fraction of the A-star lensing events (A-stars should supply ~ 1/3 of the 
events), provided that the spectroscopic observations can be initiated fast enough to 
catch the event before maximum light. The spectroscopic signature can serve as an 
additional test of the lensing (as opposed to variable-star) nature of the event. More 
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importantly, the Einstein ring radius can be measured from the temporal behavior 
of the line-shift. This technique can therefore remove part of the degeneracy in the 
derivation of Macho parameters in the ongoing experiments. 

The spectral shift can be measured with relatively smaller effort (one measurement 
per night) in a large fraction of microlensing events involving Machos in the halo of 
the LMC itself, which can then be recognized as such. If, for example, 10% of the 
events are due to LMC Machos, spectroscopic measurements of 50 A-star events with 
j3 < 0.2 for one hour per night for a week will reveal the effect in ~ 5 cases, and 
would provide a direct measurement of the LMC fraction of Macho events. The LMC 
events can be distinguished from any Galactic "foreground" events by their much 
lower proper motion (Gould 1994a). 

With a more ambitious program in which events are followed for a week with 
eight measurements per night (likely requiring an 8m telescope), the proper motion 
can also be measured for the A-star + Galactic Macho events with (3 <, 0.10. While 
this is a lot of telescope time the potential payoff is very great, particularly if the 
"reduced velocity" can be measured from parallax observations (Gould 1992, 1994b; 
Gould et al. 1994). The reduced speed is given by (-D5l ~ -^lmc)^^^^*- Hence, by 
measuring the proper motion, u;^* one would find the distance to the Macho, -Dql- 
The distance and the Einstein radius, ^*, yield the mass, M . Combining these with uo 
gives the transverse speed. In general, measurement of the reduced velocity requires 
observation from two small satellites in solar orbit. It is likely that the reduced speed 
could be obtained from only one such satellite (Gould 1994b). Gould et al. (1994) 
show that it is generally possible to measure one component of the reduced velocity 
from ground-based observations provided that the Machos are going of order the 
Earth's orbital speed of 30kms^^. This would be extremely rare for Machos in the 
halo. However, the measurement proposed by Gould et al. actually becomes much 
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more sensitive for very high-magnification events {f3 ^ 0.1) such as those for which 
proper-motion determinations can be made. For these events, we conjecture that one 
or perhaps even both components of the transverse velocity could be measured from 
the ground, even for Machos traveling at halo speeds, ~ 200 km s^^. 

Acknowledgements: We would like to thank D. Goldberg, T. Mazeh, and C 
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Figure Captions 

Figure 1 The line centroid sliift Av in a microlensed star as a function of time. The 
plots are for an Einstein-radius crossing time uj^^ = 17 days, and a projected stellar 
rotation speed tJi-ot = 100 km s '^. (The horizontal and vertical axes can be scaled 
for any uj^^ and frot)- The separate panels show the dependence of the velocity curve 
on the three parameters: dependence (a) on p, the ratio of the stellar angular radius 
to the Einstein radius; (b) on j3, the ratio of the impact parameter to the Einstein 
radius; (c) on a, the angle between the Macho velocity vector and the projected stellar 
rotation axis, a can be deduced from the shape of the velocity curve, and (3 from the 
total photometric magnification of the event. The amplitude of the velocity curve 
then yields p, or, since the stellar radius is known, the Einstein radius of the Macho. 
Figure 2 The expected relative error in p, the ratio of the stellar angular radius to 
the Einstein radius, as a function of the time span of the observations, T in days. The 
two panels are for p = 0.001 (typical for O.IMq Machos in the halo of the Galaxy) and 
for p = 0.01 (typical for Machos of this mass in the LMC halo). The various curves 
give the dependence of the error on the dimensionless impact parameter f3. Since a 
fraction j3 of lensing events have impact parameter < j3, the fraction of lensing events 
where Sp/p can be measured to a specified accuracy given a stretch of telescope time 
T can be read off the plots. It is assumed that the timescale of the event is uj~^ = 17 
days and the precision of the radial velocity measurement that can be obtained in 
1 hour at the time of maximum magnification is Sv/vr:ot = W~^. Eight hours of 
observation per day are assumed in the p = 0.001 panel and just one hour (or one 
measurement with equivalent precision) per day in the p = 0.01 panel. About 10% 
of Galactic Macho events and 20% of LMC Macho events involving A stars can be 
measured to better than 30% accuracy in 8 days of observing per event. 
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